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ABSTRACT 



Context. The compact association Cygnus OB2 is known to contain a large population of massive stars, but its total mass is currently a matter 
of debate. While recent surveys have uncovered large numbers of OB stars in the area around Cyg 0B2, detailed study of the optically brightest 
among them suggests that most are not part of the association. 

Aims. We observed an additional sample of optically faint OB star candidates, with the aim of checking if more obscured candidates are 
correspondingly more likely to be members of Cyg 0B2. 

Methods. Low resolution spectra of 9 objects allow the rejection of one foreground star and the selection of four 0-type stars, which were later 
observed at higher resolution. In a subsequent run, we observed three more stars in the classification region and three other stars in the far red. 
Results. We identify five (perhaps six) new evolved very massive stars and three main sequence 0-type stars, all of which are likely to be 
members of Cyg OB2. The new findings allow a much better definition of the upper HR diagram, suggesting an age ~ 2.5 Myr for the 
association and hinting that the 03-5 supergiants in the association are blue stragglers, either younger or following a different evolutionary 
path from other cluster members. Though the bulk of the early stars seems to belong to an (approximately) single-age population, there is ample 
evidence for the presence of somewhat older stars at the same distance. 

Conclusions. Our results suggest that, even though Cyg 0B2 is unlikely to contain as many as 100 0-type stars, it is indeed substantially more 
massive than was thought prior to recent infrared surveys. 



^ Key words, open clusters and associations: individual: Cyg 0B2 - stars: formation - stars: early-type - stars: mass function 



1. Introduction 



Among Galactic OB associations, Cyg OB2 is special in many 
respects. For a start, it is known to host a large population of 
massive stars, including a s ignificant fraction o f the earliest 
spectral types in the Galaxy (IWalborn et al.ll2002h . The optical 
extinction to Cyg OB2 is high, but not sufficiently so that it pre- 
vents spectra of its stars in the classification region being taken 
(something impossible for other very massive open clusters 
with a large population of massive stars, s uch as Westerlund 1 
dClarketal.! 12005) or the Arches Cluster dpiger et al.ll2002l) ). 
Because of this, Cyg OB2 is a very useful laboratory, since, on 
one hand, it provides a lar ge homogeneous populatio n of OB 



These properties have led to a great deal of in- 
terest in Cyg OB2, from the "classical" study of 
Johnson & Morgan (1 9541) to the c omprehensive investi- 
gation by Massev & ThompsonI (Il991 ), who identified ~ 60 



stars that can be analysed (iHerrero et al.lll999l l2002h and, on 



the other, can be used as a template to co mpare optical and in- 
frared investigations (e.g.. iHansonl 2003 ). Finally, because of 
its compactness and high stellar content, Cyg OB2 seems to 
occupy a more or less unique position somewhat intermedi- 
ate betweeii an op en cluster and a normal OB association (cf. 
Knodlsedej|2000l) . 



stars more massive than 15Mq. More recently, b ased on star 
counts in the 2MASS observations of the region, iKnodlsedei 
(I2OOOI) proposed that the number of O-typ e stars in Cyg 
OB2 was much larger Building on this result, IComeron et alJ 
(2002) preselected a large number of possible OB members 
of Cyg OB2 from their 2MASS colours and obtained low- 
resolution H- and /T-band spectroscopy of the candidates. 
Candidates that lacked molecular bands were selected as very 
likely early-type stars. Of 77 candidates so selected, 31 stars 
for which optical spectra existed were OB stars, suggesting 
that most, if not all, of the other 46 objects were also OB stars 
in Cyg OB2. 



From this list of candidates, HansonI ( 2003 ) selected those 
brightest in the optical (14 objects with B - 12 to 14), for 
which she obtained classification spec tra, find i ng tha t all of 
them were indeed OB stars. However, iHansonl (l2003h argues 
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that most of these objects are not members of Cyg OB2. For 
a start, they all lie at some distance from the previously de- 
fin ed boundaries of Cyg OB2, as most of the sources located 
by IComeron et al. (2002) do. Moreover, about half of the ob- 
jects observed are late O and early B supergiants, in dicating 
ages rather larger than the 2 Myr that HansonI (2003) derives 
for Cyg OB2 from isochrone fitting to the location of the main 
sequence. Finally, one star (A39, B2 V) appears far too bright 
for its spectral type and is almost certainly a foreground object. 

It is therefore an open question as to whether the list of can- 
didates from Comeron et al.l (l2002b really contains a high frac- 
tion of actual Cyg OB2 members. Here we investigate this issue 
with new spectra of several other fainter optical candidates. We 
also make use of the recent publication of a larg e catalogue of 
accurate spectral types for Cyg OB2 members dKiminki et al.' 



2007h which combined with our results and those of Hanson 



(120031) . allows an enormous improvement in the characterisa- 
tion of the HR diagram for the association. 



In what follows, we will use the notation ofl Comeron et al 



(120021) for stars within their list (A## for OB candidates and 
B## for emission-line star s). For other members , we w ill use 
the numbering system of Massey & ThompsonI (1991), with 
prefix M T, except for the twelve sta rs with the classical num- 
bering of Johnson & Morgaiil ( 1954 ). which are given with the 
symbol # followed by their number. 

2. Observations 



Candidate stars from the list of lComeron et alJ (120021) were ob- 
served with the 1.52-m G. D. Cassini telescope at the Loiano 
Observatory (Italy) during the nights of 2004 July 15 - 18. 
The telescope was equipped with the Bologna Faint Object 
Spectrograph and Camera (BFOSC) and an EEV camera. We 
used grism #3, which covers 3300-5800A with a resolution of 
~ 6 A. Unfortunately, on the night of July 16th, the sky was 
very poor, with some veiling, and we resorted to observing two 
stars with the lower-resolution grism #4. The night of July 17th 
we could not observe. Therefore, in total, we observed only 10 
stars, of which one, A27, had been observed before w i th bet ter 
resolution and signal-to noise ratio (SNR) by HansonI ( 2003 ). 

From these ten objects we selected five to be observed at 
higher resolution; the four which appeared to be O-type stars 
(based on the analysis presented in Section 13.1b and one that 
was probably a B-type dwarf as a check. These objects were 
observed with the 4.2-m William Herschel Telescope (WHT) in 
La Palma (Spain), equipped with the ISIS double-beam spec- 
trograph, during a service run in June 2006. The instrument 
was fitted with the R300R grating and MARCONI2 CCD in 
the red arm and the R300B grating and EEV#12 CCD in the 
blue arm. Both configurations result in a nominal dispersion of 
0.85A/pixel (the resolution element is approximately 3 pixels 
in the blue an 2 pixels in the red). 

As the selection criteria of Comeron et al. I (I2OO2') proved 
sound, we then selected some objects from their list with very 
bright K magnitudes (which should be intrinsically brightest) 
and observed them during a run on 2007 August 21-22 at the 
WHT. Three objects were observed in the blue with grating 
RI200B (nominal dispersion of ~ 0.23 A/pixel) and three oth- 



ers (whose Z? > 16 made too faint for the blue grating) were 
observed with the red arm and grating R600R in the /-band, 
wher e relatively accurate classification is also possible (e.g., 
Clark et al]|2005 ). This configuration has a nominal dispersion 
of ~ 0.5A/pixel. 

All the spectra have been reduced with the Starlink pack- 
ages ccDPACK (Draper et al. 2000) and figaro ( S hortridge e tal 



1997) and analysed using figaro and dipso (IHowarth et al 
119981) . 

3. Results 

3.1. Loiano spectra 

The Loiano spectra have rather poor SNR in the blue, but allow 
a rough classification of the stars. One of the candidates, A40, 
turns out to be a foreground G-type star The other 9 objects are 
very obviously OB stars. Their spectra are displayed in Fig.[T] 
while their 2MASS magnitudes and derived spectral types are 
listed in Table[I] 

Table 1. Infrared 2MASS photometry and derived spectral 
types for programme stars. 



Name 


U - K^) 


^s 


Spectral type 


Telescope 


All 


1.19 


6.64 


07.5 Ib-II(f) 


L, WHTl 


A12 


1.21 


5.72 


BO la 


L,WHT2 


A15 


1.14 


6.81 


07Ib(f) 


L,WHT2 


A18 


1.07 


8.35 


~08V 


WHT3 


A24 


0.97 


7.46 


06.5 III((f)) 


L,WHT1 


A25 


1.01 


7.36 


-08111 


WHT3 


A26 


0.97 


8.19 


09.5 V 


L, WHT2 


A27 


0.97 


5.75 


-BOP 


L 


A30 


0.81 


8.61 


-B2V 


L 


A31 


0.95 


7.98 


-B0.5V 


L 


A33 


0.87 


8.60 


B0.2V 


L, WHTl 


A35 


0.81 


8.47 


-BOV 


L 


A38 


0.85 


8.56 


08 V 


L, WHTl 


BIO 


1.45 


8.12 


Be 


WHT3 


B17 


1.21 


6.44 


Ofpe 


L,WHT1,WHT3 



BO la ( lHansonl l2003') 
Key for telescope configurations: 

L - Loiano Cassini Telescope, WHTl - WHT in 2006 with blue 
aim, WHT2 - WHT in 2007 with blue arm, WHT3 - WHT in 
2007 with red arm (/-band only) 



A 1 1 has He I 447 1 A =^ He n 4542A, no visible He 11 4686 A 
(at this resolution; we see it in the WHT spectrum) and 
C III 5696A strongly in emission. It is thus an -07 supergiant. 
A24 has He 1 4471 A < He 11 4542 A and He 11 4686 A in absorp- 
tion and so it is a relatively unevolved mid O-type star. A27 has 
very prominent C iii 4650A no He 11 4686 A, weak He 11 45 12A, 
strong He i lines and very weakened HjS. It should be a ~B0 su- 
pergia nt, and indeed it has been classified as BO la bv lHansoii 
(120031) , based on higher quality spectra. 

The spectrum of A30 has lower resolution and SNR than 
the rest, but extends into the red. The lack of He 11 45 1 2A makes 
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5600 5800 



Fig. 1. Useful section of the spectra obtained from Loiano. The 
SNR decreases quickly towards the blue and there are essen- 
tially no counts bluewards of ~ 4200A.The spectrum of A40, 
which turned out to be a foreground star, is not shown. 



it later than BO, while the fact that Ha is much deeper than the 
He I 6678, 7065A lines suggests that it is a mid B star (e.g., 
later than ~B2V). 

A31 has moderately strong Cm 4650A and very weak 
Hen 4686 A and Hen 45 12 A, suggesting a main sequence star 
in the BO-1 range. A33 is similar, with a slightly stronger 
C m 4650A, perhaps suggesting a higher luminosity. A35 has 
stronger Hen 4686A, but is unlikely to be much earlier, as 
Hen 4512A is weak. 

A38 has moderately strong Hen lines, but Hei 4471 A 
> He n 4542A, suggesting a late O-type star, while the lack of 
emission Unes indicates a low luminosity. Finally B17 is char- 
acterised by strong emission lines of He n 4686A and N ni, and 
may be an extreme Of supergiant or an Ofpe/WNL star. 



3.2. WHT spectra 

Figure |2] shows the spectra of the 5 objects observed in 
2006. The spectrum of B17 is very striking, with very strong 
He n 4686A and N ni emission, and a P-Cygni profile in H/3. 
All its Unes are displaced by > 200 kms ' with respect to other 




4400 4600 
WtiTcleiigtli (A) 



Fig. 2. Classification spectra of the 5 objects observed with 
ISIS on the WHT in June 2006. 



. A12 - BOIa 

/ I 1 1! ^ 1^ ^ \^ l! J 



A26 - 09. 5t 



4400 
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Fig. 3. Classification spectra of the 3 objects observed with the 
blue arm of ISIS on the WHT in August 2007. The most rele- 
vant lines are marked. 



members and show an enormous shift in radial velocity with re- 
spect to the 2004 spectrum. We classify this object as an Ofpe 
star, almost certainly a binary, and will study it in detail in a 
future paper Al 1 and A24 have He n 4542A^ He 1 4471A and 
are therefore close to 07, while the in-filling of Hen 4686A 
indicates a moderate luminosity. All has Nm 4630-4640A 
in emission, a wind featur e typical of luminous stars. B ased 
on the criteria laid out by Walborn & Fitzpatrickl ( 1990l) . we 
classify Al 1 as 07.5 Ib-II(f) and A24 as 06.5 III((f)). A38 has 
He n 4686A strongly in absorption and we classify it as 08 V, 
though i t is clo se to 08.5 V, if we use the quantitative criteria of 
Mathvl (Il988h . Finally, A33 has weak He n 4686A and 4542A, 
but no He n 4200A, and thus we classify it as BO. 2 V. The accu- 
rate classifications agree quite well with the estimates obtained 
from the low-resolution spectra in the previous section. 

Figure [3] shows the 3 classification spectra obtained in 
August 2007. The extremely prominent Si iv lines in A12 show 
it to be a luminous supergiant, while their ratio to Si ni lines 
puts it at BO, in agreement with the presence of three weak 
Hen lines. We adopt BOIa. A15 is similar to All and A24. 
He n 4686A is more clearly in emission, but the lack of wind 
S IV emission lines and weak Si iv 4089A prevent us from as- 
signing a high luminosity. We settle for 07 Ibf. Finally, though 
clearly an O-type star because of the strong He n lines, A26 still 
shows many weak Si m and O n lines and is therefore 09.5 V. 
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We also obtained red spectra of 4 objects, covering the at- 
mospheric window between 830 0A and 8900A, where spectral 
classification is possible (e.g.. iClark et al.ll2005b . Fig. |5] shows 
the spectra of three objects, B17 and two stars not observed 
in the blue, A18 and A25. The three spectra are similar and 
typical of late O-type stars. Only a few Paschen lines are vis- 
ible, but the C III 8502A line is clearly visible. This places the 
stars in the 07-09 range. The broadness of the Paschen lines 
and the small number visible indicates that A18 and A25 are 
not supergiants. We take an approximate spectral type ~08. 
Based on their K magnitudes and positions in the HR diagram 
(Section [4.31 ). A25 is likely a main-sequence object, but A18 
could be a more evolved star. 

A fourth object, BIO = MT 285, was also observed and 
its spectrum is displayed in Fig. |6l The strong asymmetric 
emission Paschen lines are typical of a Be star, with the very 
prominent Oi 8446A indicates that it is not a late-B object 
dAndrillat etaPl 19881) . The likely detection of Hei 8779A and 
strong Oi 777 4A (not shown) em ission indicates that it is 
B2 or earlier ( Andrillat et al. 19881). This i s fully consistent 
with the detection by Comeron et al. I (12002*) of He: 2.058Aim 
in emission, as this i s only seen in Be stars earlier than B3 
dClarkfe Steele! l2000h . Therefore it is likely to be a massive 
Herbig B e star in Cyg OB2. No ne of the classical Be stars ob- 
served bv lAndrillatet all d 1988b shows EWqj < -5.5, but BIO 
has EWqi =s -10 after correction for Pa 18, strongly suggest- 
ing that it is a Herbig Be object, as they tend to have stronger 
emission features. 



3.3. Model fits 

The WHT spectra, even if of moderate resolutio n, off'er a good 
chance to co mplement the study presented by iHerrero et al.l 
(Il999[ l2002h . expanding the sample that can be analysed. For 
all stars with blue WHT spectra except B17, which is unlikely 
to be a s ingle star, we determined stellar paramete rs using 
FASTwiND dSantolava-Rev et al.|[l997t IPuIs et al. l '2005*), by fit- 
ting H and He line profiles in the standard way (Herrero et al.l 
1992HRepolust et al.ll2004. The value of Voo wa s adopted from 
the spectral type, after i Kudritzki & Pulsl d2000h . We adopted a 
value of the micro turbulence ^ - lOkms"' for all objects ex- 
cept A12 (the only B-supergiant in the sample) for which we 
adopted 15 kms"'. /3 (the exponent of the velocity law) was 
adopted to be 0.8 and varied when the fit could be improved. 
Again, only A12 needed a slightly larger value of /3 (consistent 
with a slower wind acceleration), and we adopted yS = 1 .0 for 
this source (but note that lacking Ha or having low resolution, 
our data are not very sensitive to j3). We should indicate that the 
final fit to the He n 4686 A line of A12 is not satisfactory. 

Likewise, the He abundance by number relative to H plus 
He, e, is set initially to the standard value e = 0.09 for all stars, 
and varied to obtain better fits. Only A15 needed a higher value 
to fit the observed spectrum. The high value required by A15, 
e - 0.25, points to an evolved object, consistently with its low 




4460 4480 4500 

WAVELENGTH 



4520 



4540 



Fig. 4. The spectrum of A15 (07Ib; solid line) compared to 
that of Al 1 (07.5 Ib-II(f), dotted line). The strong N iii 4511 - 
14A feature and the high value of e derived from the model fits 
for A15 indicate advanced chemical evolution. The high rota- 
tional velocity and the underluminosity of this star all indicate 
anomalous evolution, perhaps due to mass transfer in a close 
binary. 



gravity and its strong N spectrum. Fig. H] shows a comparison 
of the spectra of All and A15 (both 07 supergiants), around 
4500A, where we can see the strong Nm 4511-14A feature 



in the spectrum of A15 . Note also that the projected rotational 
velocity is remarkably high (245 km s"') for an evolved object 
(which is assumed to have lost significant amounts of angular 
momentum). The values derived suggest that the evolution of 
this object has been anomalous (perhaps as a consequence of 
binary evolution) as, in addition, it appears underluminous and 
undermassive for its spectral type. 

The parameters derived are given in Table|2]and correspond 
very well to the spectral types derived in most cases. Errors 
in the stellar parameters are estimated at STf^g = +1500 K, 
6logg - +0.2 and d{logM) - ±0.3 for the low resolution 
observations and slightly lower in Teff (+1000 K) and \ogg 
(±0.15) for the higher resolution data. Absolute luminosities, 
radii and masses have been calculated fro m the A"s magn itude, 
assuming DM = 10.8 (d = 1.4 kpc), after ' HansonI d2003h . fol- 
lowing the method discussed in Section l431 In the case of A38, 
the value My - -3.7 is more than half a magnitude fainter than 
expected for the spectral type, resulting in the low derived mass 
and luminosity. 

Note that, in order to compare the absolute astrophysical 
parameters derived for thes e objects to those in previous works 
dHerrero et al. I ll 999112002 '). they must be reduced to the same 
distance, as previous works assumed the canonical DM = 11.2. 

Projected rotational velocities could not be determined for 
the stars observed in 2006 due to the low resolution; the instru- 
mental profile dominates the line broadening for the metals (H 
and He lines are broadened by the Stark profile). However, this 
provides an upper limit, as instrumental broadening dominates 
in all our objects, allowing us to ascertain that they all rotate 
with V sin i < 160 km s ' . 
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Name 


Spectral type 


Mv 


TeS 


log,? 


R 


M 


I'oo 


Vrot 


Mass 


log(L/Lo) 








(K) 




(Ro) 


(Mq) 


(kms-') 


(km s-') 


(Mq) 




All 


07.5 Ib-II(f) 


-5.8 


36000 


3.6 


15.9 


2.2 X 10-« 


1900 


< 160 


38.9 


5.6 


A12 


BO la 


-6.7 


27000 


3.0 


30.2 


3.5 X 10-* 


1350 


80 


34.2 


5.6 


A15 


071bf 


-5.7 


35000 


3.2 


15.6 


3.2 X 10-« 


2100 


245 


19.0 


5.5 


A24 


06.5 lll((f)) 


-5.0 


37500 


3.6 


10.7 


1.7 X lO-'' 


2600 


< 160 


18.1 


5.3 


A26 


09.5V 


-4.2 


35000 


3.9 


7.7 


4.1 X lO^** 


1300 


90 


17.6 


4.9 


A33 


B0.2V 


-3.6 


31000 


4.0 


6.6 


2.0 X lO^** 


1000 


< 160 


16.6 


4.6 


A38 


08V 


-3.7 


36000 


4.0 


6.0 


4.9 X 10^** 


1900 


< 160 


13.8 


4.7 



All models have been calculated using /3 = 0.8, Vtmb = 10 km s ' and e = 0.09, except for A12, which, being a B-type supergiant, required 
different wind parameters (fi = 1.0, V(„,.i, = 15 km s"') and A15, which required e = 0.25. 



cm e5Da(+pi6) 




Wavelength (A) 



Fig. 5. /-band spectra of three targets observed with the WHT 
in August 2007. A25 does not show clear evidence for any 
Paschen line beyond Pa 1 3 and has weak C m 8502 A. It is hence 
most likely an 08-9 V star. A18 is clearly more luminous and 
probably earlier. Note the important radial velocity shifts in the 
lines of B17, this time in the opposite sense to those in the blue 
spectrum in Fig. |2] 




0] li'l'iS P16 850a P15 B545 Fi4 aSSB P13 966 



ffavelengtli (A) 

Fig. 6. /-band spectrum of B 10. The data available identify it as 
an early (< B2) Be star, but do not allow us to decide whether it 
is a classical Be star or a PMS Herbig Be object. The enormous 
strength of 1 8446A, though, is unusual for a classical Be star 
and points to the second option. 



4. Discussion 

4.1. Completeness 



Our results confirm the enormous success of Comeron et al 



(120021) at identifying reddened OB stars. Only one of the can- 
didates turns out to be an interloper. The important point. 



logt=6.4 (rot) 




Fig. 7. Semi-observational HR diagram for Cyg OB2, based on 
published spectral types and 2MASS JHK^ photometry (see 
text for details). Continuous lines are non-rotating isoch rones 



The 



for log? : 

dashed line is the log? - 6.4 is ochrone in the high-rotation 
models dMevnet & Maedejl2003h . 



however, is estimating whether these objects are members of 
Cyg OB2. The line of sight in this direction runs parallel to the 
Local Arm, and populations at different distances may lie pro- 
jected together. While it is extremely unlikely that early O-type 
might be found far away fro m massive cltisters or associations, 
except for a few runways (cf . de Wit et al.l2005 ; van den Bergh 



20041) . less massive stars will certainly be found if one looks 



through a Galactic Arm. 

In this s ense, our sample appears rather different from that 
of HansonI (2003), who observed only candidates which were 
bright in B, most of which turned out to be foreground B-type 
stars. Even though this is not surprising for the August 2007 
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sample, which was selected on the basis of bright K magni- 
tudes, it is more striking for the first (Loiano) sample, which 
simply consists of objects somewhat fainter in B than those ob- 
served by Hanson (2003). This suggests that, even though ex- 
tinction is clearly variable across the face of Cyg OB2, on aver- 
age, there is a range of extinctions where we can find members, 
and this translates into a range of magnitudes (according to 
spectral type). This ran ge has been estimated as 4 ^ Ay ^ 7 by 
previous authors (e.g., iMassev & Thompsonlll99ll) . and may 
extend to somewhat higher values when stars from the list of 
IComeron et al. I (l2002h are added. 

Because of this, we suggest that there cannot be many 
more O-type memb ers amongst the candidates given by 
Comeron et al. (l2002h . though there must be some (for exam- 
ple, perhaps A17, see below). Very few of them are likely to be 
intrinsically bright (and so very massive, evolved) members of 
the cluster. Based on their { J -K^) colours and magnitudes, 
only A4 and A8 might be sufficiently bright (intrinsically) to 
be obscured O-type giants or supergiants. 

In addition, our data reveal three new evolved O-type stars 
(All, A15 and A24), which help define the main-sequence 
turn-off of the association. All is of particular interest, as it 
lies very close to the Blue Hypergiant (BHG) candidate #12, 
in what likely is the most obscured part of the association. We 
classify it 07.5Ib-II(f), as it almost looks evolved enough to 
be a supergi ant, and its analysis in deed shows that it is a very 
massive star ^ Comeron et al" 1 (12002 ') suggest it may be the coun- 
terpart to th e X-ray source IE 202 3043-H4103.9. Another can- 
didate from lComeron et aP ( 2002 ). A17, lies very close to it. It 
has very similar IR colours, but is two magnitudes fainter in K. 
This is most likely a late-O/early-B main-sequence member. 

4.2. Clustering 

Fig |8] shows a 2MASS image of the central region of 
Cy g OB2, containin g the two cluster-like groupings identified 
by iBica et al. I (12003*) and the area around #12. The two clus- 
ters are prominent against the background. The field shown in 
Fig.[8]is ~ 12' X 12', corresponding to ~ 5 pc at 1.4 kpc. The 
separation between Cluster 1 and Cluster 2 (^ 6') is equivalent 
to 2.5 pc, and so smaller than the radius of relatively massive 
clusters in the Perseus Arm, such as h Per or NGC 663. 

The nine brightest stars in Cluster 2 have (J - K)^ = 
0.63 ±0.03 (standard deviation) and eight of them have {J-K)^ 
between 0.59 and 0.63. This uniformity in reddening repre- 
sents strong confirmation of their associ ation, also c lear i n 
the colour-magnitude diagram shown by iBica et alj (l2003b . 



Cluster 1 , on the other hand, does not seem to have a uniform 
reddening, but we do not think that this is a strong argument 
against its reality, in this region of patchy obscuration. 

In addition to these two groups, and outside the field cov- 
ered by Fig. [8] there is another o bvious region of stella r over- 



density - present in the data of iKiminki et alJ (l2007h - sur 



rounding star #4 (07 III). It comprises MT213 (BO V), MT215 
(B2V), MT216 (B1.5 V) and MT221 (B2V). MT187 (Bl V), 
MT227 (09 V), MT241 (B2 V) and MT258 (08 V) He within 
3'. MT187 and MT221 are significantly (> 0.1 mag) more red- 



dened than the others, but the other seven have (J - K)s - 
0.49 + 0.03, again strongly hinting at a real physical associa- 
tion. 

The presence of all these small groups over a large area 
suggests that star formation has proceeded in small bouts in 
this region, perhaps over an extended period of time. In spite 
of this, the bulk of the population occupies positions in the HR 
diagram incompatible with a very long period of star forma- 
tion. Subclustering is seen i n the largest Galactic s tar form- 
ing re gions, such as W49A (iHomeier & Alvesll2005h or W51 
dNand a-Kumar et a l.ll2004l) . Study of large stellar complexes in 
M5 1 ( Bastian et alj|2005b shows that the age spread within dif- 
ferent clusters is ^ 10 Myr. 



4.3. HR diagram and ages 



As discussed bv lHanson ( 2003b . the main sequence in Cyg OB2 
ext ends clearly down to 06V. Star #22, classified 04 111(f) 
by iMassev & ThompsonI (1199 11) has been shown to be a 
close double containing an 03 If* supergiant and an 06 V 
sta r (Walborn et al. 2002 ) MT 516, classified as 05.5V((f)) 
by IMassev & ThompsonI ( 199 lb. was found to have a rather 
low gravity by HerreroetalJ ( T999I) . Indeed, the fact that 
He II 4542A is somewhat stronger than He ii 4686A shows that 
this object is rather far away from the ZAMS. Therefore this 
star is likely better classified as 05.5 III, joining #8C and the 
faint component of #8A as an object still on the main sequence, 
but already showing some signs of evolution. Therefore the age 
of the association would seem to be set by the fact that stars 
more massive than 06 V are already somewhat evolved, while 
06 V stars are not. 

However, within a classical theory of stellar evolution, it is 
difficult to see how this fits with the presence of 03 supergiants. 
In order to address this question and also exploit the potential 
of Cyg OB2 as a laboratory, we have constructed an HR dia- 
gram utilising the wealth of new spectral type determinations 
in this region. We have used the 2MASS JHKs ma gnitu des for 
all objects and their spectral types (from iKiminki et al.ll2007i 



Hansonll2003l or this work) in order to place them in a semi- 



observational HR diagram. We have followed the procedure 
used by, for instance, IMassev et al.l (Il995h . but taking infrared 
rather than optical magnitudes. We have resorted to JHKs mag- 
nitudes partly because many stars of interest lack good f/-band 
photometry, but also because this allows a test of the usefulness 
of infrared data to study obscured massive clusters. 

From the spectral types derived, we have taken a 
Teff an d bolometric correc tion BC, using the calibra- 
tion of Martins et al. ' ('2005') for O-type stars and that of 
ISumphrevs & McElrov. (.1984) for B-type stars (the two cal- 
ibrations agree quite well around BO; however, the possible 
existence of an artificial jump between BO and B 1 has been 
noted by previous authors). We also take intrinsic (V - K )o and 
(J - K)() colours from the calibration of IWegnerl (Il994l) . With 
the observed (J-Ks), we derive E(J-Ks). As the reddeni ng to 
the as sociation is known to be very close to standard tHansonl 
2003b . we simply calculate Ak^ = 0.67£'(7 - Ks). 
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Fig. 8. A 2MASS A's image of the central region of Cyg OB2. All catalogued stars more massive than BO V are marked. MT573 
(B3 1) is obviously a background object, while the connection of MT556 (Bl I) and MT601 (BO lab) to the association is unclear. 
They may repres ent an older (~ 8 Myr) population at about the same distance. The two cluster-like groupings identified by 



Bica et al.l (120031) stand out in the image. Object 1 is centred on MT425 and contains #22, while Object 2 is centred on the three 



components of #8. 



\ ye then calcula te Kq^ Ks- Ak, - DM, using DM = 10.8 
from lHansonI (l2003h . and by adding (V-K)o and the BC, arrive 
at a semi-observational Mboi. Fig.|7]plots Mboi against the Teff 
derived from the spectral classification. Superposed on it, are 
Geneva isochrones without rotation for log t - 6.2 (1.5 Myr), 
log t - 6.3 (2 Myr) and log t - 6.4 (2.5 Myr), as well as the ro- 
tating isochrone for log f = 6.4. We tried to fit the data using the 
higher distance modulus (DM = 1 1 .3) obtain ed by averaging 
spectroscopic distances (IKiminki et al.ll2007h . but this left all 
the stars well above the ZAMS. The data used for this diagram 
are hsted in Table 3. 



The diagram shows several notable features. There is a 
very well tra ced main seque nce extending to the 06V stars. 
As noted by iHansonl (120031) . A37 (05 V) is earlier than any 
other MS stars, but seems to fit the sequence well. Around 
log Teff ~ 4.3, there lie a number of B stars well above the 
main sequence, unlikely to be connected with the rest of the 



population. Objects like MT642 (B 1 III), A23 (B0.7 lb) or A34 
(BO. 7 lb) are very probably not members of Cyg OB2. 

There are two evolutionary sequences that seem to turn off 
the main sequence. The one below the isochrones is formed by 
MT138 (08I), A32 (09.5IV), A41 (09.7II) , A29 (09.7Iab), 
A36 (BO lb), and MT601 (BO lab). Though this sequence may 
provide a decent fit to the log t - 6.7 (5 Myr) isochrone, the 
random distribution in luminosity class suggests that this is 
not a real evolutionary sequence, but simply the projection of 
a number of luminous stars situated at slightly different dis- 
tances. It is worthwhile mentioning, though, that many of these 
objects have distance moduli comparable to that of the main 
Cyg OB2 association. 

In contrast, the sizable population of evolved stars lying 
around the logf = 6.4 isochrone seems to form a much more 
coherent group. Moving along the isochrone, we have A24 
(06.5III) and #4 (07III), #8B (07II-III), All (07.5Ib-II), 
A20 (08II), #10 (=MT632, 09.5Ia), A12 (BO la) and A27 
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(BO la). The excellent progression in luminosity class with 
spectral type strongly supports the hypothesis that these objects 
are really following the isochrone. Only two stars with accu- 
rate spectral types do not fit this evolutionary sequence: one is 
MT771 (07V), which appears as bright as #4. This is easily 
explained by the fact that it is a double-li ned spect roscopic bi- 
nary with two similar components ( Kimi nki et al ]r2007). The 
other one is A15 (07Ib), which is ~ 0.6 mag fainter than ex- 
pected. As mentioned, the analysis of its spectrum reveals very 
high He and N abundances, and a very low mass for its spec- 
tral type. Therefore this is indeed a peculiar object, perhaps the 
product of mass transfer in a close binary, as #5 and B17 are 
also likely to be. 

The distribution of main sequence stars, the main-sequence 
turn-off and the sequence of evolved stars strongly supports an 
age ~ 2.5 Myr for the bulk of Cyg OB2. According to the 
caUbration of IMartins et al.l (l2005b . this implies that stars up 
to ~ 35 Mq are still close to the ZAMS, while more massive 
stars are akeady more evolved. However, it is obvious that the 
brightest stars in the association fall well above the adopted 
isochrone. 

Given the extent of Cyg OB2, the possibility of a spread in 
ages cannot be e xcluded and may e ven seem logical. Indeed, 
in a recent paper, i Drew et aL ( 2008 ) have shown an important 
concentration of A-type stars to the South of the bulk of the 
O-type stars. In order to be at the same distance as the O-star 
association, these A-type stars must be part of a 5 - 7 Myr pop- 
ulation. 

Does the presence of 03 supergiants indicate the existence 
of an even younger population? Certainly this possibility can- 
not be excluded, but it is worth taking in consideration two 
points: 

- If there is an age difference, we would expect to find some 
sort of spatial segregation between the older and younger 
population, but this is not evident in the data. The area 
shown in Fig. |8] contains most of the earliest objects, but 
also Al 1 (07.5 Ib-II) and #12. The moderately evolved #8B 
(07 II-III) falls just in the middle of Cluster 2, which con- 
tains three of the early objects. 

- Stars more massive than the 06V objects still in the 
main sequence appear as either earlier-type supergiants 
(Or stars) or intermediate luminosity 06-7 stars. In other 
words, if there is a younger population, all its members ap- 
pear as or stars just now. 

At the estimated age, it is unlikely that any stars might have 
undergone supernova, and indee d no supernova remnant is seen 
in the area (iPasquali et al ] l2002h . All the stars in Fig. |7] occupy 
positions in the HR diagram compatible with being still in the 
hydrogen core-burning phase. The BHG candidate #12 may be 
past this phase and there are 5 Wolf-Rayet star s in the region 



that h ave been proposed as possible members (Pasqual i et al 



locked in an LBV phase, like # 12) is highly suggestive of the 
idea that not all very massive stars evolve in the same way. 

There is, however, ample evidence suggesting that star for- 
mation has been going on for quite some time in a large area 
around the recognizable core of Cygnus OB2. Indeed many of 
the evolved massive stars that are unlikely to belong to the cur- 
rent generation of massive stars lie at approximately the same 
distance and could belong to an older (~ 7 Myr) generation, 
associated with the young A-type stars detected by 'Drew et al.l 
(2008 ). In the massive association 30 Dor. Walb orn & Bla deSj 
(Il997h also find a population of OB stars somewhat older (4 - 
6 Myr) than those in the main cluster R136 scattered across 
the entire complex. Likewise, Mokiemet al. (2007) find ages 
of 7.0 + 1.0 and 3.0 + 1.0 Myr for the associations LH9 and 
LHIO, in the giant H ii region Nil in the LMC. Their data are 
consistent with LH9 having triggered star formation in LHIO. 



5. Conclusions 



20021). The WC stai's WR 144 and WR 146 may actually be 
the descendants of the most massive stars in the association 



(iPasquali et al.ll2002l) . The fact that some of the most massive 
stars appear as 03 If supergiants, while others are moving to- 
wards the red part of the HR diagram (or seem already to be 



Though the candidate sample of IComeron et al.l (l2002h con- 
tains a high fr action of likely non-members, as discussed by 
HansonI (2003), it has also allowed the detection of a number 
of obscured O stars and very luminous BO la supergiants very 
likely to be members of Cyg OB2. 

When these objects are included in the HR diagram, it be- 
comes clear that there is a sequence of moderately evolved stars 
detaching from the main sequence exactly at the position where 
we stop seeing luminosity class V objects, i.e., around 06 V. 
These two facts combined support an age of ~ 2.5 Myr for the 
bulk of the association. 

The HR diagram presented in Figure [T] contains the largest 
number of Cyg OB2 members ever displayed in such a dia- 
gram. It contains ~ 50 stars that may have started their lives 
as main-sequence O-type stars and only a few of these are un- 
likely to be members. Unless a population of extremely ob- 
scured O-type stars is lying at fainter magnitudes than probed 
by 2MASS, the total number of O-type stars in the association 
is almost certain to be in the 60-70 range. 

The number of stars that have already left the main se- 
quence and lie above the 06 V members that define the turn- 
off is more securely determined. If the main association is ba- 
sically co-eval, these represent the subset of stars that were 
originally more massive than 35 Mq. Counting #12, which is 
not shown in Fig. Ill bec ause of its claimed spectral variabil- 
ity dKiminki et al.ll2007h . there are 21 such stars. The evolved 
interacting binaries #5 and B17 (not in Fig. |7]i should be 
counted too (perhaps doubly). The resulting number is cer- 
tainly only a lower limit. Apart from possible unrecognised 
close doubles and binaries, at least two of the Wolf-Rayet stars 
in the area are likely to be descendants of very massive stars 
(IPasquali et al.ll2002l) . Also. IComeron & Pasqualil (l2007b iden- 
tify BD -i-53°3654 as a likely runaway 04 If member of the as- 
sociation. Therefore, we have identified a population of at least 
25 stars that were originally more massive than 35 Mq. Given 
the uncertainties - in particular the very high binary fraction 
dKiminki et alj|2007i) - we refrain from trying to derive a total 
mass for the association by assuming an IMF. 
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The brightest members, with spectral types in the 03-05 
range, may technically be considered blue stragglers. Though 
a real age difference cannot be ruled out, it does not seem to be 
borne out by the spatial distribution of stars, perhaps suggesting 
that we are seeing stars of similar mass evolve in very different 
ways. 

Three luminous supergiants (#10 09.5Ia, A12 BO la and 
A27 BO la) seem to follow the 2.5 Myr isochrone and so appear 
to be the descendants of stars more massive than ~ 40 Mq. This 
is in agreeme nt with an initial mass estimate of 48 Mq for #10 
( Herrero et al.ll2002l) . which may have to be slightly reduced if 
the lower DM - 10.8 is adopted. These objects will probably 
soon reach the LBV instability, which #12 is perhaps already 
encountering. A large population of 09-B1 la supergiants de- 
scended from stars with M, * 35 Mq is found in the older 
(~ 4.5 Myr) cluster Westerlun d 1 together with a number of 
LBVs and Yellow Hypergiants ( Clark et al.ll2005h . 

In summary, even if Cyg OB2 falls short of the proposed 
100 O-type stars by a factor of ~ 2, its nuclear region still rep- 
resents one of the most massive collections of early-type stars 
known in the Galaxy and its relatively low reddening cements 
its value as a laboratory for the study of their properties. 
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